Evolution of condensate fraction during rapid lattice ramps 
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Combining experiments and numerical simulations, we investigate the redistribution of quasi- 
momentum in a gas of atoms trapped in an optical lattice when the lattice depth is rapidly reduced. 
We find that interactions lead to significant momentum redistribution on millisecond timescales, 
thereby invalidating previous assumptions regarding adiabaticity. We show that this phenomenon 
is driven by the presence of low-momentum particle-hole excitations in an interacting system. Our 
results invalidate bandmapping as an equilibrium probe in interacting gases. 



Introduction — Optical lattice experiments are making 
important contributions to our knowledge of equilibrium 
and non-equilibrium properties of quantum many-body 
systems [l| . A prime example is that of bosons in opti- 
cal lattices Q , which are described by the Bose Hubbard 
Hamiltonian. There are many of the open questions in 
this system related to transport and dynamics Here 
we explore one such issue: how does the condensate frac- 
tion evolve during a lattice ramp? 

In addition to its relevance for developing paradigms of 
nonequilibrium physics in cold gases, our study is of prac- 
tical value: lattice ramps are routinely used in a probe 
known as bandmapping [4|-l6j . The bandmapping protocol 
involves three steps. First, the lattice depth is reduced 
over a carefully chosen timescale r, with the intention 
of mapping quasi-momentum onto momentum. Second, 
the gas is allowed to ballistically expand for a short time, 
mapping momentum onto position Q. Third, the gas is 
imaged. These images have been interpreted as the initial 
quasi-momentum distribution. This protocol has been 
used to measure condensate fraction [8(, map Brillouin 
zones 0], determine temperatures and probe phase 
transitions 1 l|-tl3| . Here we show that bandmapping is 
an inherently non-equilibrium process. Understanding 
the timescale of evolution of quasi-momentum is there- 
fore crucial to interpreting the data. We find that in 
interacting systems the dynamics are much faster than 
previously believed, leading to systematic errors in mea- 
surements of quantities such as condensate fraction. 

The key step in bandmapping is turning off the 
lattice slowly enough such that the quasi-momentum 
states adiabatically evolve into momentum states, yet 
fast enough to leave the occupations of different quasi- 
momentum states unchanged. For harmonically trapped 
non-interacting particles, there are three relevant en- 
ergy scales: the bandgap Eb g , the tunneling t, and the 
quantum energy of the harmonic confining potential hv 
(h = 2irh is Planck's constant). Provided that the lattice 
ramp time r is long compared to the inverse bandgap 
(r ^> h/Ebg), the quasi-momentum states will adiabati- 
cally evolve into momentum states (lol | . To avoid motion 
of atoms in the trap, the lattice must be turned off quickly 
compared to the trap period (r <C 1/v)- Thus, h/Ef, g 



and 1/v set the adiabatic and diabatic timescales respec 



tively. Aside from its influence on Eb g and v (14|, the 
hopping rate t has little direct impact on bandmapping. 
In standard experiments in the tight-binding regime, the 
natural separation of timescales between h/E^g < 0.1 ms 
and 1/v > 10 ms makes it straightforward to satisfy the 
adiabatic condition; a ramp time r ~ 1 ms is usually 
employed. 

Complications arise for an interacting gas. Interactions 
lead to a coherent redistribution of quasi-momentum oc- 
cupations [Ulllfij]- Furthermore, collisions scatter atoms 
between different quasi-momcnta while conserving total 
quasi-momentum. The on-site interaction energy be- 
tween two atoms U determines the relevant dynamical 
timescale h/U; in most experiments, h/U < 1 ms. We 
show that the additional criterion, r <C h/U, disrupts the 
separation of timescales that makes bandmapping suc- 
cessful in noninteracting systems. 

Using a combination of experiment and numerical sim- 
ulations, we investigate the impact of interactions on 
bandmapping for atoms confined in a lattice in the 
strongly correlated regime. We quantify the redistri- 
bution of quasi-momentum during lattice ramps for a 
Bose-Einstein condensate of atoms in a 3D cubic opti- 
cal lattice. The fraction of atoms in the condensate is 
determined after linearly ramping from lattice depth Vi 
(with IOEr < Vi < 14Er, spanning the supcrfluid and 
Mott- insulator regimes) to a fixed final depth Vf = 4JEr. 
Here En = (h/X) 2 /2m is the recoil energy, A is the laser 
wavelength, and m is the atomic mass. The final depth 
is chosen so that the atoms remain in the single-band 
Bose-Hubbard limit, thereby simplifying comparison to 
numerical simulations. Ramps terminating at Vf = 
produce similar results. Our results apply to interacting 
gases in general, including fermionic systems and mix- 
tures. 

Experimental Method. — Our experimental setup is de- 
scribed in detail in Ref. [Io| . In summary, we create a 
condensate composed of 87 Rb atoms in the \F — 1, mp = 
— 1) state in a harmonic trap with (geometric) mean trap 
frequency D = 35.78(6) Hz. We cool the gas until the 
condensate fraction exceeds 80%. 

We superimpose a cubic optical lattice with a d = 
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A/2 = 406 nm lattice spacing on the atoms by slowly 
turning on three pairs of retro-reflected laser beams. 
The laser intensity determines the potential depth V. 
Through Kapitza-Dirac diffraction, we calibrate V to 
within 1%, but drift in the calibration results in a 6% 
systematic uncertainty. The Gaussian envelope of the 
lattice beams adds to the harmonic confinement, and the 
overall (geometric) mean trap frequency with the lattice 

(120 ± 10) nm is 
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■, where w 



the measured 1/e radius of the lattice laser beams. 

Ten milliseconds after loading the lattice, we linearly 
ramp V from V to Vf = 4 En in time r. The lattice and 
trapping potentials are then removed in 10 ns and 0.2 ms, 
respectively, and the column density is imaged after the 
gas expands for 20 ms. We extend the dynamic range 
of our measurement by imaging only a controlled frac- 
tion of atoms that are transferred to the F = 2 hyperfine 
state. The number of condensate atoms N c is measured 
using multimodal fits to "low optical density" images for 
which only a small number of atoms are transferred. We 
supplement these with "high optical density" images, for 
which all are transferred and imaged. In these images, 
the broad non-condensate component is resolvable, but 
the condensate peaks are saturated (see Fig. [T]). The 
number of non-condensate atoms is determined by fitting 
the broad background with the condensed peaks masked 
and extrapolating the non-condensate component into 
the masked regions. The total number of particles N var- 
ied from (103 ± 5) x 10 3 to (72 ± 2) x 10 3 for V = 10 E R 
to 14 Er, and the condensate fraction ranged from 0.3 
to 0.05. 

Experimental Results. — For each V, we measure the 
post-ramp condensate fraction as a function of r. A typ- 
ical data set is shown in Fig. [TJfor = 10 En. The data 
points follow an exponential, as illustrated by the red 
curve. The fitted time constants r rc i are shown in Fig. [2] 
as a function of the initial condensate fraction (bottom 
axis) along with the corresponding Vi (top axis). The re- 
laxation time is relatively weakly dependent on Vi, only 
changing by a factor of two as the superfluid-Mott insu- 
lator transition is crossed, and throughout is consistent 
with the simple empirical rule r rc \ ocl/U. 

The insets to Fig. [T] shows the high optical density 
images for short (a) and moderate (b) ramps. Even for 
ramps as short as 1ms, the quasi-momcntum distribution 
changes dramatically, as atoms are transferred to low mo- 
mentum. Most bandmapping experiments use ~ 1 ms for 
r, and therefore do not measure the initial condensate 
fraction. Shorter ramps arc not a solution to this prob- 
lem, as they lead to significant non-adiabatic transfer of 
atoms to higher energy states [17|. 
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FIG. 1: Condensate fraction measured after bandmapping 
from Vi = 10 Er. The insets show high optical density im- 
ages where the background is resolved, but the Bragg peaks 
are saturated. The images are shown in false color, with red 
(blue) indicating regions of high (low) column density. The 
field of view is 813 x 813 /im. (a) r = 10 ns; (b) r = 1 ms. 
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where at and a\ are bosonic annihilation and creation 
operators at lattice site i, and /i, = /i — V e -x.{i), where /i 
is the chemical pot ential and V ex (i) is the external po- 
tential at site i [l9(. The first sum in Eq. 1 is over all 
nearest neighbor sites. We calculate U using the exact 
Wannier functions in the lowest band j^, and extract 
the tunneling amplitudes t from Mathieu characteristics. 

We calculate dynamics using a time dependent 
Gutzwiller ansatz |18| . which approximates the wave- 
function by * = (g),- J2 m c ™ ( t )l m )» wnere \m)i is the 
m-particle Fock state on site i, and the coefficients cf (t) 
are space (i) and time (t) dependent. (Note the typo- 
graphic distinction between the tunneling t, time t, and 
ramp time r.) This approximation leads to a simplistic 
quasi-momentum distribution, dividing atoms into zero 
momentum (k = 0) condensed and k ^ non-condensed 
states. The total number of condensed atoms N c is given 

by N c =Eil(ai)| 2 , where (oj) = £ m yfr+l<$ ~ W 

As described in Ref. 
ihdtip = for \P yields a set of coupled differential 
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Theoretical Modeling. — We use the 3D Bosc Hubbard 



where the mean-fields are fa = £{j)( a j)/6j m which the 
sum is restricted to nearest neighbors of site i. 
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FIG. 2: Relaxation time (r re ;) for the condensate fraction 
for ramps from Vi to Vf = 4Er for variable ramp times r. 
The range of Vi spans the superfluid (SF) and Mott insulator 
(MI) regimes (demarcated by the vertical dotted line). The 
experimental data is bounded by zero temperature Gutzwillcr 
mean-field simulations using two different initial states (see 
text). The solid black line shows the relaxation time assum- 
ing initial state 2, while the dashed line shows the relaxation 
time assuming initial state 1. The error bars represent the 
uncertainty in the relaxation time from a fit to data such as 
that shown in Fig. [TJ 



As in the experiment, the potential depth is ramped 
V{t) =Vi + (V f - Vi)(t/T), where Vi and V f = 4 E R are 
the initial and final lattice depths, and r is the ramp time. 
The Hubbard parameters t and U are time-dependent 
because of this ramp. In the simulations we include a 
spherically symmetric external harmonic trapping poten- 
tial V ex , matched to the (lattice-depth-dependent) exper- 
imental value v. 

We make direct comparison with the experimental data 
by studying TV ~ 75, 000 87 Rb atoms on a 55 x 55 x 55 
lattice with lattice spacing 406 nm. For our initial state 
we use a local density approximation obtained by solv- 
ing the homogenous, single-site problem. To account for 
the overestimation of the condensate fraction in mean- 
field theory, we use two different initial states to model 
the data. Initial state 1 is the mean-field ground state 
obtained by using the physical lattice depth Vi in the sim- 
ulation. This state has a larger condensate fraction than 
the experimental initial state. Initial state 2 is obtained 
by finding the lattice depth at which the condensate frac- 
tion predicted by the theory matches the measured con- 
densate fraction at Vi. 

Time evolution, from either of these initial states, is 
calculated using a split-step approach with sequential site 
updates. We ramp the lattice down from Vi to AEr in 
a time r ranging from to 1.5 ms, and calculate the 
condensate fraction at the end of each ramp. We then 
fit the resulting data to an exponential curve to extract 
a characteristic relaxation time T re /. 

Comparison of Theory and Experiment - In Fig. [2] we 
compare two theoretical curves with the experimental 
data. 



The bottom line is obtained from initial state 1, and 
the top line is obtained from initial state 2. Initial state 
1 yields a higher condensate fraction compared to the 
experimental system at any given initial lattice depth. 
It relaxes to equilibrium faster, hence providing a lower 
bound on the characteristic relaxation time r re ;. Initial 
state 2 treats the atoms as if they were in a deeper lattice 
than the physical system, and therefore leads to slower 
dynamics. It therefore provides an upper bound on the 
characteristic relaxation time. The dashed line indicates 
the lattice depth for which a shell of unit-filling Mott- 
insulator emerges. Throughout, the data points lie be- 
tween these two theoretical bounds. 

For Vi < 13 En both theoretical protocols yield sim- 
ilar results. Here the entire system is superfluid, and 
mean-field theory is accurate. Throughout this regime 
TYci ~ 0.5ms. As the Mott transition is approached, the 
relaxation time increases by a factor of ~2, indicative of 
slower dynamics in the insulating state. The simulation 
using initial state 2 (top curve) captures this physics, 
showing a significant increase in relaxation time; initial 
state 2 contains a Mott-insulator shell. The growth of 
the relaxation rate from initial state 1 is more gradual, 
as the Mott-insulator transition occurs for larger values 
of Vi as compared to the experiment. While the slower 
timescales for dynamics in the Mott-insulator state are 
quite intuitive, the restoration of phase coherence follow- 
ing a rapid quench is not fully understood 22. 23|. 

Our simulations used zero temperature initial states — 
finite temperature would modify the connection between 
condensate fraction and lattice depth, effectively raising 
the relaxation rate obtained from initial state 1. 

Understanding the timescales. — Bandmapping 
timescales are too short for any significant transport in 
the lattice to occur 
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and thus the observed physics 
is purely local. One gains insight by considering the 
homogenous case V ex = 0. Linearizing Eq.J51 about the 
homogenous Gutzwillcr stationary state |24j, we find 
the excitation spectrum uj(k) for the lattice gas in the 
shallow and deep lattice limits (Fig. [3j . 

In the shallow lattice, the only relevant excitations are 
linearly dispersing phonons. Excitations that change the 
overall phase of the wave-function cost no energy (c.f 
Fig |3] (left)). On short timescales, the overall density re- 
mains constant and phonons are not excited. For deeper 
lattices, a particle-hole branch with a quadratic disper- 
sion, and a gap A is present. 

The superfluid state in a deep lattice is "quantum de- 
pleted": the fraction of condensed (k = 0) atoms ap- 
proaches near the insulating transition. Upon lower- 
ing the lattice depth quasi-momentum states coherently 
evolve from high momentum to low momentum. Full 
phase coherence is restored when these excited atoms 
subsequently decay into particle-hole pairs of compara- 
ble energy but low momentum, while changing the phase 
of the local wave-function. The timescalc for this pro- 
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FIG. 3: Typical Gutzwiller excitation spectra for deep lat- 
tices. (Left): Quantum depleted superfluid near the Mott- 
insulator transition (V = 13-Bh). Two modes are present: 
a gapless phonon mode with a linear Bogoliubov (dashed) 
dispersion at low k, and a gapped particle-hole mode with 
quadratic dispersion. The gap A ~ U, sets the diabatic- 
ity timescale for bandmapping. (Right): Mott insulator at 
(i = U/2 and V = 28-Er. The lowest energy excitation has a 
gap at k = of order U/2, that sets the diabaticity time. 



cess is set by 1/A, the timescale to excite a particle- hole 
excitation. 

We produce an analytic expression for A by truncating 
the Fock basis to at most 2 particles per site and diago- 
nalizing the resulting 3x3 Hamiltonian. The eigenvalues 
become particularly simple when \i w U/2, and we find 
that A = \{U + y/A8z 2 (a)H 2 + U 2 ), where z is the co- 
ordination number (z = 6 for a cubic lattice). At t — > 
and unity filling, the gap A — > U/2 ~ 1.5ms in deep lat- 
tices (Fig. [3]). Therefore, even in the superfluid phase, 
interactions foil bandmapping because it is not possible 
to ramp off the lattice quickly compared to h/A and 
slowly compared with h/Eb g . 

Summary. — We find that previous assumptions that 
dynamics are frozen during bandmapping are incorrect 
for interacting systems. Rather, we observe that con- 
siderable momentum redistribution occurs for typical 
bandmapping times, driven by the presence of particle- 
hole excitations in the quantum depleted superfluid. 
These excitations render bandmapping unreliable for 
measuring quasi-momentum distributions in interacting 
systems. This problem could be circumvented by turning 
off interactions prior to bandmapping using a Feshbach 
resonance. 

We remark that much of the physics of these ramps is 
captured by time-dependent Gutzwiller mean-field the- 
ory. There is, however, room for quantitative improve- 
ment: for ramps from deeper lattices, different assump- 
tions lead to a five-fold variation in the relaxation time 
T re ;. To accurately model the experimental findings, as 
well as answer subtler questions, such as how correla- 
tions propagate across the gas, one must resort to more 
sophisticated techniques. [251 ] . 

Extracting information about the underlying many 
body physics from quench experiments is a vibrant area 
of research 



lattice ramp is an indirect measure of particle-hole exci- 
tations. Probes such as bandmapping are therefore in- 
dispensiblc to developing our understanding of strongly 
correlated systems. 
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whose importance transcends ultra- 
cold atoms. Here we have shown that the timescales gov- 
erning the evolution of the condensate fraction during a 
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